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ABSTRACT 

We test and reject the claim of Segal et al. (1993) that the correlation of 
redshifts and flux densities in a complete sample of IRAS galaxies favors a quadratic 
redshift-distance relation over the linear Hubble law. This is done, in effect, by treating 
the entire galaxy luminosity function as derived from the 60-^m 1.2 Jy IRAS redshift 
survey of Fisher et al. (1995) as a distance indicator; equivalently, we compare the flux 
density distribution of galaxies as a function of redshift with predictions under different 
redshift-distance cosmologies, under the assumption of a universal luminosity function. 
This method does not assume a uniform distribution of galaxies in space. We find that 
this test has rather weak discriminatory power, as argued by Petrosian (1993), and the 
differences between models are not as stark as one might expect a priori. Even so, 
we find that the Hubble law is indeed more strongly supported by the analysis than 
is the quadratic redshift-distance relation. We identify a bias in the the Segal et al. 
determination of the luminosity function, which could lead one to mistakenly favor 
the quadratic redshift-distance law. We also present several complementary analyses 
of the density field of the sample; the galaxy density field is found to be close to 
homogeneous on large scales if the Hubble law is assumed, while this is not the case 
with the quadratic redshift-distance relation. 

Subject headings: 

1. Introduction 

While the precise value of the Hubble constant Hq is still the subject of controversy (cf. 
Jacoby et al. 1992), the linearity of the redshift-distance relation is generally not questioned by 
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most astronomers. Since its original announcement by Hubble (1929), observational evidence has 
mounted steadily in its favor; for recent discussions, see Lauer & Postman (1992), Peebles (1993), 
Riess, Press, Sz Kirshner (1996), and Perlmutter et al. (1996). However, Segal and collaborators 
have persistently argued for a quadratic redshift-distance relation wherein z oc r p , where z is the 
redshift, r is the distance, and p = 2 (Segal et al. 1993, hereafter SNWZ; Segal & Nicoll 1995, 
and references therein). This relation is predicted for the low-redshift regime in the Chronometric 
Cosmology developed by Segal (1976). 

The greatest difficulty both in establishing the Hubble Law and in measuring Hq is 
determining the distances to astronomical objects; this is often accomplished through the use of 
standard candles. The approach of Segal and collaborators in effect is to use the entire luminosity 
function of galaxies as a standard candle. If one assumes that the luminosity function of galaxies 
is independent of position and local density, at least after averaging over solid angle in a large 
redshift survey, then the comparison of the derived luminosity function as a function of redshift in 
principle could be a test of the assumed redshift-distance relation. Equivalently, the distribution 
of apparent magnitudes as a function of redshift in principle contains information about the 
redshift-distance relation (but see Petrosian 1993 and below). SNWZ choose to examine low-order 
statistics from the observed distribution of apparent magnitudes m and redshifts z in a given 
sample. In particular, they calculate the scatter in m, and the slope of the regression of magnitude 
on log redshift, comparing the expected values (under the assumption of various values of p) to 
what is actually observed. The distribution of apparent magnitudes or flux densities at a given 
redshift is independent of the density distribution of galaxies, and makes no assumption about the 
large-scale homogeneity of the universe. SNWZ have carried out this analysis using the 1.936 Jy 
redshift survey of galaxies observed with the Infrared Astronomical Satellite (IRAS) of Strauss 
et al. (1990; 1992). They found that p = 1 (the Hubble law) was strongly rejected by their analysis, 
with p = 2 (the Lundmark (1925) law) reproducing the observed magnitude-redshift correlations. 

Despite the importance of Segal's claims over the years, there has been very little response 
in the literature to this work. Soneira (1979) calculated the quantity (z\m) for the galaxies in 
the Reference Catalogue of Bright Galaxies (de Vaucouleurs & de Vaucouleurs 1964) and found 
that p = 1 was greatly favored by the data. However, this analysis requires the assumption of a 
homogeneous galaxy distribution, unlike the analysis of SNWZ (cf., Nicoll & Segal 1982). A similar 
analysis has been carried out by Shanks, as quoted in Peebles (1993, §5); see also Choloniewski 
(1995). A more general discussion of the consequences of the quadratic redshift-distance law is 
given in Salpeter & Hoffman (1986). Finally, recent comparisons of observational data with the 
predictions of Chronometric Cosmology at high redshift (which we do not address in this paper) 
can be found in Segal & Nicoll (1986), Wright (1987; cf., Segal 1987), Cohen et al. (1988; cf., Segal 
1990), Efron & Petrosian (1992), and Segal & Nicoll (1996). 

Petrosian (1993) argues that given flux densities and redshifts for a complete sample, it is 
impossible to separate cosmological effects (including evolution) from the luminosity function; one 
must make additional assumptions about underlying distributions. He therefore argues that one 
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cannot use the observed correlations between flux density and redshift as a test of cosmologies. 
The results of this paper are in accord with this; we find that without additional assumptions 
about the underlying density field, the tests carried out by Segal are only very weak discriminants 
of cosmological models. 

In the present paper, we carry out an analysis similar to that of SNWZ, using data from a 
redshift survey of IRAS galaxies. We discuss the derivation of the universal luminosity function 
in §2. In §3.1, we follow SNWZ in using the distribution of redshift and flux density as a function 
of redshift as a test of cosmological model. In §3.2, we allow ourselves the assumption that the 
galaxy distribution is homogeneous on large scales, giving us a variety of further cosmological 
tests. We conclude in §4. 



2. Derivation of the Luminosity Function 

The luminosity function <&(-£/) of galaxies is the distribution of galaxies as a function of 
luminosity: <&(L) dL is the mean number of galaxies per unit volume, with luminosity between L 
and L + dL. We will convert from redshift to distance by writing cz = h p r p , and arbitrarily adopt 
values of h p = 100 km s~ x Mpc~ p . With these conventions, the luminosity L of a galaxy is given 
by L = Air(cz/h p ) 2 ^ p uf, where / is the observed flux density at frequency v. 

Following the notation of Yahil et al. (1991) we will find it useful to define the cumulative 
luminosity function ^(L): 



*(L)= /<&(L')dL' , (1) 

implying that 



L 



d^ 

The cumulative luminosity function is closely related to the selection function (f>(z), which is the 
fraction of the luminosity function which enters the sample at a given redshift. For a sample that 
is flux-density-limited to / > f m in, we have 

> y[L min (z)\ 

where 

Lmin = 4:TT(cz/h p ) 2/p uf min (4) 

is the minimum luminosity detectable at redshift z, and L m i n (z s ) is a self-imposed lower limit on 
luminosity, corresponding to that of a galaxy at the flux density limit placed at cz s = 500 km s _1 . 

Methods for deriving luminosity functions are reviewed by Binggeli, Sandage, & Tammann 
(1988) and by Strauss & Willick (1995), and can be classified as parametric or non-parametric. 
The former assume some implicit functional form for the luminosity function, characterized by a 
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number of parameters which are adjusted to provide the best fit to the observational data. The 
latter make no assumptions as to the form of the luminosity function, and are constructed directly 
from the data. However, both approaches assume that the luminosity function is universal; that is, 
independent of position or local density. This assumption of universality is of course a necessary 
one if we wish to use the observed luminosity distribution for cosmological tests (although, as 
Petrosian 1993 argues, and as we show below, it is not sufficient). 

Here we will derive the luminosity function by maximum-likelihood estimation, following 
Nicoll & Segal (1978; 1980), Sandage, Tammann, & Yahil (1979), Nicoll & Segal (1983), Efstathiou, 
Ellis, & Peterson (1988), and Yahil et al. (1991). If the luminosity function is universal, then the 
joint probability density that we have a galaxy of luminosity Li with redshift Zi is a separable 
function of these two variables: 

F(Li,Zi)dLdz = $(L i )p(z i )dLdz , (5) 

where p(z) is the local density of galaxies at redshift z. The conditional probability density 
T{Li\zi) that the galaxy have luminosity Lj, given its redshift Zi, is then given by the joint 
probability, divided by the integral of the joint probability over all possible luminosities in the 
sample at that redshift, given the flux density limit. That is, 

T{Li\zi) = = (6 ) 

iZLfe) H L > z i) dL J $(L)dL 

Note that the density field p(z) has dropped out of this equation. The quantity L m i n is defined in 
Eq. (|j); the p-dependence enters the analysis when calculating the luminosity of each galaxy from 
its redshift and flux density. 

We then maximize the likelihood of observing the entire sample: 

C = Y[F(Li\zi) . (7) 

i 

For computational convenience, we instead minimize the negative logarithm of the likelihood, 
where 

-ln£ = -^ln^(Li|zi) . (8) 

i 

An attempt was made to include p as one of the parameters with respect to which C was 
maximized, but this proved impossible; C turns out to be a monotonically increasing function of p. 

The maximum-likelihood approach does not depend on the assumption that galaxies are 
distributed uniformly, since p(z) dropped out of Eq. (0). However, this requires that the mean 
spatial density of galaxies in the sample be calculated by other means. There is a variety of such 
density estimators available (Davis & Huchra 1982); the one chosen here is 

n 1 = y[L min (z s )] = V- 1 J2^-s . (9) 
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where the sum is over all galaxies in the sample volume V, and 4>{zi) is the value of the selection 
function at the redshift of the i th galaxy, as in Eq. (||) . 

The data set consists of 5313 pairs of 60 /xm IRAS flux densities and redshifts above a flux 
density limit of 1.2 Jy, taken from Strauss et al. (1992) and Fisher et al. (1995). The recession 
velocities are corrected for the motion of the Sun with respect to the barycenter of the Local Group 
following Yahil, Tammann, & Sandage (1977). Of these galaxies, 4218 have redshifts in the range 
500 < cz < 12000 km s -1 ; only these were used in fitting the luminosity function since the sample 
may be incomplete at much higher redshifts (Fisher et al. 1992), and at low redshifts the local 
motions are strongly affected by peculiar velocities that dominate the local Hubble expansion. 



2.1. Parameterized derivation 

Various parameterizations of the luminosity function are discussed by Strauss (1989) and 
Saunders et al. (1990). We follow Yahil et al. (1991), in parameterizing the cumulative luminosity 
function by 

m{L) = c(^X a (l + ^-Y , (10) 



^ J \ L 
so that the differential luminosity function is given by 



dL \L L* + L 

and the minimization is performed with respect to the two dimensionless parameters a and (3 and 
the characteristic luminosity L*, using the routine mrqmin from Press et al. (1992). The optimal 
values of the parameters for p = 1, 2, 3 are tabulated in Table H; the case p = 3 will be included in 
some of the tests we present for methodological perspective. 

For graphical purposes, it is often convenient to work with the distribution per volume per log 
luminosity, given by 

^ = - h/* n = T^mL) ■ (12) 
d(logio-L) logioe 

We shall use this form in making plots of the luminosity function. 



2.2. Non-parameterized derivation 

One drawback of assuming a parameterized form of the luminosity function is that it 
constrains the luminosity function to have a certain functional form, which may not be general 
enough to accurately reflect the actual luminosity function. A non-parametric approach to 
determining the luminosity function does not suffer from this limitation. 
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We model the differential luminosity function as being piecewise-constant over n bins evenly 
spaced in log 10 L from the minimum to maximum luminosities seen in the sample, and then treat 
the value of the function on these bins as the n parameters with respect to which the likelihood 
(defined exactly as before) is to be maximized. This avoids any implicit assumptions of the 
functional form that the luminosity function should have, but also sacrifices any requirements as 
to continuity and smoothness that physical intuition suggests should be satisfied. This approach 
was first suggested by Nicoll & Segal (1980; 1983), and was reinvented by Efstathiou et al. (1988). 
Rather than attempting to search for a minimum in some high n-dimensional parameter space, the 
technique is to converge iteratively to the optimal step values. The interested reader is referred to 
the Appendix or Efstathiou et al. (1988) for details. 

There are several drawbacks to the non-parametric luminosity function method. The most 
obvious is that the resultant luminosity function is discontinuous; this can be overcome by 
interpolating linearly between the centers of the bins; the details are set forth in the Appendix. We 
in fact use this interpolating method in what follows below. More serious is the effect illustrated 
below; a strong sensitivity of the derived luminosity function to bin size. 

The ROBUST method of SNWZ calculates the luminosity function in two steps. For a given 
sample with maximum redshift z max , there is a luminosity L m i n (z max ) = Air (cz max /h p ) 2 ^ p vf m i n 
above which a galaxy can be found anywhere in the volume of the survey, so &(L) in this 
luminosity range is simply proportional to the number of galaxies at each value of L. For lower 
luminosities, however, the volume in which a galaxy could be found is an increasing function of 
the luminosity, requiring further calculation to determine the luminosity function. Nicoll & Segal 
(1983) and SNWZ argue that the number of bins nu used in the lower-luminosity range alone 
should be set the same when comparing different cosmologies, in order to avoid giving any one 
cosmology extra degrees of freedom. Clearly, nu < n by definition. 

With this in mind, Figs. |l] and ^ show the resultant luminosity functions derived under the 
assumption of p = 1 and p = 2, respectively. We plot the parameterized luminosity function 



(Eq. 11, using the parameters of Table |T]) as a solid curve, and superimpose the luminosity 
functions derived using the nonparametric method for several different binnings: the bin size used 
by SNWZ (in logL), nu = 10 as recommended by SNWZ, and nu = 250. We plot error bars only 
for the non-parametric luminosity function that most closely approaches the parametric one. The 
bin sizes and number of bins for these and other luminosity functions are tabulated in Table § 
below. We also plot the luminosity functions of SNWZ for comparison (normalized using Eq. |9|). 
Different bin sizes result in appreciable differences in the faint end of the luminosity function; 
in particular, the fewer the number of bins, the more the faint end of the luminosity function is 
attenuated. 



3 Because the range of luminosities below Lmi n (zmax) is dependent on the sample and the assumed value of z„ 
our nu — 10 case has a different bin size than that of SNWZ. 
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At the bottom of these plots we have placed histograms indicating the distribution of 
the luminosities from which our luminosity functions were computed; note that the luminosity 
distribution of the galaxies in the sample narrows significantly with increasing p. 

For p = 1, note that the luminosity function we find using the SNWZ binning is close to the 
luminosity function found by SNWZ themselves. However, this bin size is large enough to strongly 
bias the luminosity function at the faint end; the luminosity function does not become stable 
until the number of bins approaches n = 20 or more (corresponding to nu = 15). Thereafter, the 
luminosity function found by the non-parametric method, which is in principle free to assume 
any shape at all, is in good agreement with the parametric luminosity function. For p = 2, the 
sensitivity to binning seems to be much less severe, and the luminosity functions we find for all 
binnings are in very good agreement. In this case, we were able to reproduce the SNWZ luminosity 
function only with a very small number of bins nu = 5, corresponding to n = 10 (not shown). 

The data set used by SNWZ consisted of the brighter half of the present one, flux density- 
limited at 1.936 Jy at 60 um. This makes little difference; we find that the luminosity function 
from different subsets in flux density is quite robust (Koranyi 1993), as long as nu is large enough. 

The biasing in the luminosity function due to binning is more severe for p = 1 than for 
p = 2, and therefore statistics for different values of p using these luminosity functions can give 
misleading results if the binning used is overly coarse. Indeed, we now show that with the bin size 
recommended by SNWZ, p = 2 is indeed a better fit to the data than is p = 1, while with finer 
binning, p = 1 is preferred. 

3. Comparison of Different Cosmologies 

Having derived a luminosity function for the sample under the assumption that the density 
and luminosity distributions of the sample are separable, we are now in a position to test the 
relative merits of the Hubble and Lundmark laws. We can test for robustness of the luminosity 
function, and also for self-consistency in the predictive powers of the luminosity function under 
the assumption of various power law cosmologies. These comparisons separate naturally into those 
that do and do not depend on the assumption that the galaxy distribution is uniform on large 
scales. We start with the latter, following SNWZ. 

3.1. Density-independent comparisons 

SNWZ argue strongly that it is important to distinguish between cosmological tests that 
involve quantities which depend on the cosmology itself (such as luminosity, or the luminosity 
function), and those that are pure observables (such as flux density, or equivalently magnitude, 
and redshift). Moreover, they develop statistical tests that do not depend on the assumption that 
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the distribution of galaxies is uniform in space. The statistics on which they put the greatest 
weight are s, the standard deviation of apparent magnitudes in a sample, and /?, the slope of the 
regression of the apparent magnitudes on log redshift (the notation is that of SNWZ). Before 
calculating these statistics, it is useful to examine Fig. [3|, which is the observed relation between 
apparent magnitudes (defined here, following SNWZ, as m = 60 — 2.51og 10 /) and log 10 c,2. The 
correlation between these two quantities is not very strong. The line shown is the regression line 
of m on log 10 cz for the subsample of galaxies with 2000 < cz < 20,000 km s _1 ; we have found 
that these statistics are very sensitive to the exact lower-redshift cutoff at redshifts below 2000 
km s _1 . Table § gives the observed values of s and (3. SNWZ calculate the expected values of 
these statistics using Monte-Carlo simulations. Here we will do so analytically. 

Given a model for the luminosity function and a value of p, one can calculate the probability 
distribution function of flux density / of galaxy i, given its redshift Z{. The distribution function 
for a single galaxy is given by the luminosity function, normalized appropriately, following Eq. (|6|): 

where (f>(z) is the selection function, L = A'Ki'(czi/hp) 2 / p f , and L m in was defined above in Eq. (Q). 
From this distribution function, we can easily calculate the moments of the apparent magnitudes, 
and therefore the scatter s. Similarly, we can calculate the expectation value of m for each 
value of z, from which /3 follows directly. Note that both these statistics depend on the observed 
distribution of redshifts, and cannot be considered independent of this. 

Table ^ shows the predicted results for p = 1,2, and 3 for various binnings for the subsample 
of galaxies with 2000 < cz < 20,000 km s _1 . It is clear that the results are quite sensitive to the 
details of the binning. Let us start by concentrating on the results using the binning of SNWZ 
(the first row for each value of p). Both p = 2 and p = 3 do a much better job of predicting the 
observed apparent magnitude scatter s than does p = 1. We interpret this to be due to the fact 
that the p = 1 luminosity function is much more biased at this coarse binning than that of p = 2. 
We find similar results at n\\ = 10, as recommended by SNWZ. None of the models do a good job 
of reproducing the slope of regression, 0, with this binning, with p = 1 overpredicting the observed 
value of (3 by as much as p = 2 underpredicts it. It is clear from Fig. |3], however, that a linear fit 
to the magnitude-log redshift scatter diagram is a rather poor way to model the data, and that 
the results may be quite sensitive to a small number of outlying points. 

With finer binning (jiu = 25), and with the parametrized luminosity function, the predictions 
for s are essentially independent of p. As Petrosian (1993) first argued, and as we conclude below, 
without additional assumptions, it is very difficult to distinguish cosmologies from redshift and 
magnitude data alone. However, note that the p = 1 does predict the correct value of (3 with this 
larger number of bins (the parametrized model does less well), while p = 2 and p = 3 fail quite 
badly. We describe the x 2 column in this table below. 

We can ask more of the data than simply these low-order statistics. Indeed, Eq. (13) gives a 



- 9 - 



prediction for the flux density distribution of galaxies as a function of redshift; we can compare 
this directly with what is observed. That is, the sum of these distribution functions over some 
subsample of a redshift survey can be compared with the observed flux density distribution as 
an a posteriori test of the luminosity function (Sandage et al. 1979; Yahil et al. 1991; Strauss 
&; Willick 1995), or, in the present application, of the cosmology assumed. In other words, the 
predicted distribution of galaxies of a given flux density / is 



where the summation is over all galaxies in a particular (sub)sample. The predicted and actual 
values in each bin can be compared with the x 2 statistic, yielding a measure of how similar the 
two distributions in fact are. 

Fig. |i] shows the results for the full sample of galaxies (500 < cz < 20,000 km s _1 ), as well 
as for a variety of redshift subsamples. The observed distribution in flux density is indicated by 
the dots, and the predictions for p = 1 and p = 2, using the luminosity functions with nu = 25 
are indicated by the solid and dashed curves, respectively. As expected for a flux density-limited 
sample, the distributions peak strongly towards the flux density limit of 1.2 Jy, with a tail of 
higher flux density observations; since this must be the case regardless of cosmology, we are not 
surprised that there are no gross morphological differences between the flux density distributions 
for the competing cosmologies. Indeed, the difference between the predicted curves for p = 1 and 
p = 2 is very small. Following Yahil et al. (1991), we calculate the x 2 statistic of the difference 
between the predicted and observed curves, using Poisson error bars and summing only over bins 
with five or more galaxies. The results are tabulated in Table ||[ The fit with p = 1 is acceptable 
in all bins, p = 2 fares rather worse; although it is acceptable in several of the redshift bins, the 
fit for the full sample is unacceptable. For p = 3 (not shown in the figure) the fit in most redshift 
bins is unacceptable. 

The quantity v in Table || is the number of bins of flux density in which the comparison is 
done. One might argue that this number should be reduced by the number of free parameters 
in the luminosity function (25, in these cases!). We make two points here: first, we have found 
qualitatively very similar results to those presented here when we use the parameterized luminosity 
function of Eq. (pT|), which uses only three parameters (Table |l|). Second, note that the fit to the 
luminosity function is done for galaxies with redshifts between 500 and 12,000 km s _1 ; for p = 1, 
this fit remains good for galaxies between 12,000 and 20,000 km s _1 , while for p = 2 and p = 3, 
the fit is unacceptable in this range. 

Fig. [| repeats this exercise using the SNWZ binning. The fits are not nearly as good as before, 
especially at the lowest luminosities. This is quantified in the x 2 statistics tabulated in Table f|. 
However, although no value of p is acceptable, p = 2 is much preferred over p = 1, especially for 
the full sample. We saw in Figs. || and || that nu = 10 corresponds to a bin size which gives a 
strongly biased luminosity function for p = 1, but that the bias was much less severe for p = 2. 




(14) 
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We explore the effect of binning further in Table |2[ The x 2 column gives the x 2 values 
(for v = 19) for the flux density comparison for a variety of binnings, for the sample with 
2000 < cz < 20,000 km s _1 . By this statistic, p = 1 is strongly ruled out with nu = 10 or with 
the SNWZ binning (doing much more poorly than even p = 3), while with either nu = 25 or with 
the parametrized luminosity function, p = 1 gives acceptable results. 

We believe that this is the origin of the claims of SNWZ. With proper determination of the 
luminosity function, the distribution function of flux densities given the redshifts can be predicted 
almost as well with p = 2 as with p = 1; it is very insensitive to cosmology, as was pointed out 
originally by Petrosian (1993). However, the luminosity function is biased if one chooses too low 
a value of nu, and this effect is more severe for p = 1 than for p = 2, causing one to erroneously 
conclude that p = 2 is preferred by the data. 

Our conclusion from this discussion, mirroring that of Petrosian (1993), is that any statistic 
derived from the distribution of flux densities and redshifts will be able to match the data roughly 
equally well for p = 1 or p = 2, without making further assumptions about the distributions. 

One such assumption we could make is that the universe approaches homogeneity on the 
largest scales (the Cosmological Principle; cf., Peebles 1993). With such an assumption, we can 
derive further statistics that do allow a sharp distinction between different values of p. 



3.2. Density-dependent comparisons 

Under the homogeneous approximation, one can predict the distribution of galaxies with 
redshift, given the luminosity function and a value of p. In a shell of thickness Az at a redshift z, 
one expects there to be 

/ \ 3/p-l 

N(z) = -n x <j>(z) ^ Az (15) 

P \ h pJ 

galaxies, where u> is the solid angle covered by the sample (11.06 ster in our case; cf. Strauss 
et al. 1990). These predictions for various p are shown in the upper panel of Fig. [6|, along with 
the observed distribution with redshift (we use nu = 25 throughout this section). The lower 
panel shows the ratio of the observed to predicted counts (which is the fractional density of the 
shell relative to the mean over the sample). The p = 1 curve agrees closely with the observed 
distribution, although the agreement is not perfect. The effects of the Local Supercluster at 
cz ~ 1500 km s _1 and of the Great Attractor and Pisces-Perseus regions at cz ~ 4500 km s _1 are 
visible as overdensities, even when averaged over almost the entire sky. At very large redshifts, the 
p = 1 curve slightly overestimates the observed distribution, probably because of incompleteness 
of the sample at high redshifts (Fisher et al. 1992). However, if either the p = 2 or p = 3 model 
were correct, one would have to argue that as one looks to higher redshifts in the universe, the 
inhomogeneities grow. The galaxy density plotted in the lower panel would need to be a strong 
function of redshift: low nearby, rising rapidly to a maximum at 5000 - 10,000 km s~\ and then 
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dropping by a factor of two thereafter. This would violate the Cosmological Principle, and indeed, 
if such massive structures were common in the universe, we would see them reflected in the angular 
correlation function of faint galaxies (e.g. Maddox et al. 1990). 

Indeed, one can calculate the distribution of densities on shells as a function of redshift, 
without any calculation of the luminosity function at all (Saunders et al. 1990). In § 2 above, we 
calculated the distribution function of luminosities conditioned on the redshifts T{Li\zi). Here we 
calculate the distribution of redshifts conditioned on the luminosities T{zi\Li), which is given by: 



fdz%F(L i: z) ~ A^^dzzVP-ip(z) 



where 

K ( L, \P/ 2 



-uiii.i-. i - P ( . „ ) • (17) 

When we conditioned on redshift (Eq. |^), the density distribution dropped out of the expression, 
while here, the luminosity function drops out. We can now maximize the likelihood with respect 
to the density field defined at a series of steps, exactly as we did with the luminosity function^, 
using the iterative technique described by Efstathiou et al. (1988). The results are shown as points 
in the lower panel of Fig. |(| The solid points are for p = 1, the open circles for p = 2, and the 
stars for p = 3. Error bars are given only for p = 1 to keep the figure from getting overly crowded; 
the error bars are in fact quite insensitive to the value of p. The agreement between this density 
field and that given by the curves, which depends on the calculation of the luminosity function, is 
striking; this shows us that both the luminosity function and density field calculations are robust. 

A final approach is suggested by Soneira (1979) and Choloniewski (1995); cf., Nicoll & Segal 
(1982). The expectation value of logz as a function of flux density follows directly from Eq. (|l^): 

nngy | f v _ f dz% log zF(L,z) _ fdz log zz^p(zML) 

{0gZU) Jdz^T(L,z) JdzzV P -^p(z)<l>(L) ' U8j 

where L = AT:(cz/hp) 2 ^ p uf . Unlike the expressions for the conditional probabilities, neither p nor 
$ drops out of the expression. Fig. ^ shows this statistic in bins of log flux density for the full 
IRAS sample (solid points, with errors in the mean shown) and also for the Northern and Southern 
Galactic hemispheres of the sample separately (open circles and stars, respectively), as a test of 
the robustness of the statistic to density inhomogeneities. The averaging is done for galaxies in 
the redshift range 500 km s -1 to 20,000 km s _1 . The p = 1 prediction assuming p(z) = 1 is given 
as the light solid line, while p = 2 is shown with the light dashed line. These curves are not pure 
power laws, because of the upper limit on redshift we imposed. The p = 1 line is not a perfect fit 
to the data; one is presumably seeing the residual effect of density inhomogeneities, as is made 



4 We will not, however, apply the additional complication of interpolating the density field between bins as we do 
in the Appendix for the luminosity function. 
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clear by the differences between the Northern and Southern hemispheres, especially at large flux 
densities. However, p = 2 does much more poorly, especially at large flux densities. This is not 
unexpected: as Fig. || showed, the homogeneity assumption is a very poor one for the p = 2 model. 
We can include this effect by carrying out the integration of Eq. ([0]) , using the density field found 
non-parametrically in Fig. ||; the results are shown as the heavy solid and dashed line for p = 1 
and p = 2, respectively. The p = 1 line is now in excellent agreement with the data, while the 
p = 2 line approaches the data points somewhat more closely, but is still far from a good fit to 
the data. It is not clear whether this represents a breakdown of the universal luminosity function 
assumption which went into Eq. (18), for the case of p = 2. 



4. Conclusions 

The density-dependent methods of comparing the cosmologies favor p = 1 quite unambiguously. 
However, one reaches this conclusion only with the auxiliary assumption that the density field 
of galaxies approaches homogeneity on large scales. In the case of the density-independent 
comparisons of the distribution of flux density, we have recourse to quantitative x 2 testing of 
goodness-of-fit, and we find that the standard p = 1 cosmology yields better agreement than does 
the p = 2 case. But perhaps most surprising is the extent to which these methods do a poor job of 
discriminating between different values of p, implying that the results of any analysis predicated 
on their use should be interpreted with great caution. These results are in accord with the analytic 
arguments of Petrosian (1993): without additional assumptions, questions of cosmology and the 
appropriate luminosity function cannot be decoupled from flux density and redshift data, and that 
one cannot determine both simultaneously 

A further methodological pitfall which affects SNWZ is the sensitivity of the non-parametric 
luminosity function to bin size. SNWZ used very coarse binning and therefore were working with 
luminosity functions inaccurate at the faint end. Since this effect is less pronounced for p = 2 than 
for p = 1, the luminosity function that SNWZ employed for p = 2 was less in error than the one 
they employed for p = 1. We found that with equally coarse binning, this effect causes p = 2 to 
be favored, both by the statistics used by SNWZ and by our own comparisons of the flux density 
distribution. 

Therefore from redshift survey data alone, one can conclude that p = 1 is preferred only 
weakly over p = 2 if we only allow ourselves the assumption of a universal luminosity function. 
When one makes the additional assumption that the distribution of galaxies approaches isotropy 
on large scales, the case for p = 1 becomes much stronger; indeed, if p = 2, one would need to 
argue that the density field of galaxies beyond 10,000 km s _1 drops dramatically and steadily with 
redshift. When this is combined with the observed linearity of the redshift-distance diagram using 
measurements of extragalactic standard candles (cf., Mould et al. 1991; Lauer & Postman 1992; 
Hamuy et al. 1995, 1996; Perlmutter et al. 1996; Riess et al. 1996), the evidence for the Hubble 
law becomes overwhelming. 
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A. Interpolation of the piecewise-constant luminosity function 

The piecewise-constant luminosity function as defined by Efstathiou et al. (1988) has one 
drawback; the corresponding selection function (defined from Eqs. || and |l]) shows a scalloping 
effect (cf., Strauss & Koranyi 1994). This is illustrated in Fig. [8|, which shows the ratio of 
the selection function calculated from the piecewise-constant luminosity function (dashed line) 
calculated with bins of AlogL = 0.15, to the analytic selection function (from Eq. [jl]). These 
were calculated for p = 1. The selection function of the piecewise-continuous luminosity function 
has a discontinuity in slope at the edge of each bin (as indeed it must; the slope is proportional to 
the luminosity function itself). The effect is minimized here by using small bins, but nevertheless 
it would be best to eliminate this altogether. We do so by generalizing the method of Efstathiou 
et al. (1988) to define the luminosity function not as constant in a series of bins, but rather as 
a series of line segments connecting bins. The formalism of Efstathiou et al. can be carried over 
exactly, if we replace Equations 2.8-2.11 of that paper with the following. Define the luminosity 
function to be a series of line segments connecting a series of points (L^, $fc), k = 1, . . . , ttQ; thus 

n 

<S>(L) = Y / ®kW(L,k) . (Al) 

k=l 

Note that in the present paper, we refer to the luminosity function as <£, while Efstathiou et al. 
refer to it as <p. The logarithm of the likelihood is then given by 



In C = In 



N n 
i=l k=l 



In 



N n 

X! $k[H(L min (zi), fc) - H(L max (zi), k)] 
i=l fe=l 



(A2) 



where N is the total number of galaxies in the sample. Note that Eq. (2.9) of Efstathiou et al. is 
erroneously missing the summation over k in the first term on the right-hand side. Here, L m i n 
is defined in Eq. @, and L max (zi) = \ix(cz /h p ) 2 l p v f max if one derives the luminosity function 
for samples with an upper limit on flux density. If there is no upper limit on flux density, 
H(L max (zi),k) = 0. 

The expressions for W and H are more complicated than those in Efstathiou et al. (1988), 



5 Unlike Nicoll & Segal (1983), we do not explicitly split our calculation of $ into those regions above and below 
Lmin(zmax). There is therefore no reference to mi in this Appendix. 
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but they just involve linear interpolation: 

L-L k . i 

W(L, k) 



Lk—Lk-i 



Lk+i—Lk 





Lfc-l < L < Lfc and k ^ 1 

Lk < L < Lk+i and k ^ n 
otherwise. 



(A3) 



H(L,k) 



( 2(^+1 ~~ ^fc-i) 

I ifc+l — ifc + {Lk + L — 2Lk~l)j; 

1 (L fc+1 -L) 2 

2 L k+1 —L k 

I 



L fe -L 



£ < ^fe-i and & 7^ 1 
Lfc_i < L < Lk and A; 7^ 1 

Lk < L < Lk+\ and k ^ n 
otherwise. 



(A4) 



Note that here, we define L n+ i = L n . Finally, the integral constraint we use is: 



Y.®k{L k 



k-lj 



(A5) 



One can normalize after this using the quantity n% defined in Eq. (^) (as indeed we have done in 
Fig. |8|). The solid line in Fig. |8| shows the ratio of selection function for the resulting continuous 
luminosity function to that for the analytic luminosity function; the scalloping effect has gone 
away. Other than the scalloping effect, the interpolated selection function, and that using the 
piece- wise continuous method are in good agreement, implying that this interpolation technique 
has very little effect on the derived luminosity function (although both differ at the 10% level from 
the analytic luminosity function; this is a residual effect of the finite binning, as described in §2). 
However, it does have a non-negligible effect on the derived selection function, and we have used 
it in all calculations requiring a nonparametric luminosity function in this paper. 
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Table 1. Parameters of the parameterized luminosity function. 



p 


a 




log 10 (L„) 
L Q 


m 

(10- 6 ( km s" 1 )- 3 ^) 


1 

2 
3 


0.49 ± 0.07 
-0.02 ± 0.14 
-0.81 ±0.30 


1.81±0.12 
2.47±0.14 
3.07±0.24 


9.681°;!° 
7.981°:!° 
7.161°;!° 


0.058 
4.357 
25.25 



Table 2. Statistics of SNWZ for the sample with 2000 < cz < 20, 000 km s -1 . 



Model 






A log L c 


3 d 


P a 


x 2f 




Real data 








0.628 


0.610 






p=\ 


9 


13 


0.365 


0.686 


0.700 


377.7 


Bin size used by SNWZ 


p=l 


10 


14 


0.302 


0.661 


0.677 


203.0 




p=\ 


25 


37 


0.109 


0.599 


0.591 


10.7 




p=l 








0.587 


0.536 


5.1 


Parameterized Luminosity Function 


p = 2 


8 


18 


0.187 


0.616 


0.520 


115.7 


Bin size used by SNWZ 


p = 2 


10 


21 


0.137 


0.604 


0.489 


73.0 




p = 2 


25 


53 


0.052 


0.592 


0.452 


40.4 




p = 2 








0.588 


0.444 


38.9 


Parameterized Luminosity Function 


p = 3 


7 


22 


0.131 


0.613 


0.369 


103.9 


Bin size used by SNWZ 


p = 3 


10 


29 


0.089 


0.608 


0.349 


76.7 




p = 3 


25 


71 


0.034 


0.604 


0.332 


59.4 




p = 3 








0.597 


0.328 


71.8 


Parameterized Luminosity Function 



a The number of bins in the luminosity function below L m in{z max ) . 
b The total number of bins in the luminosity function. 
c The size of the bins in log L. 

d The standard deviation of the apparent magnitudes. 
e The slope of the apparent magnitude regression on logc2. 

f The \ 2 of the predicted flux density distribution of the sample to that observed. 
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Table 3. \ 2 analysis for flux density distribution of galaxies, nu = 25. 



cz 

km s" 1 


V 


# of Galaxies 


x 2 


p = \ 

X> 


x 2 


p=2 
X> 


x 2 


p = 3 

X> 


500-2000 


19 


674 


23.1 


1.21 


22.0 


1.16 


39.8 


2.09 


2000-4000 


17 


901 


14.2 


0.84 


20.5 


1.21 


23.8 


1.40 


4000-6000 


15 


1068 


13.8 


0.92 


17.0 


1.14 


19.9 


1.32 


6000-8000 


13 


697 


18.5 


1.42 


26.4 


2.03 


32.1 


2.47 


8000-12,000 


11 


878 


5.0 


0.46 


8.6 


0.78 


17.6 


1.60 


12,000-20,000 


10 


639 


9.8 


0.98 


16.6 


1.66 


24.9 


2.49 


Entire Sample 


19 


4857 


12.2 


0.64 


45.7 


2.40 


70.4 


3.70 



Table 4. \ 2 analysis for flux density distribution of galaxies, using SNWZ binning. 



cz 

km s^ 1 


V 


# of Galaxies 


x 2 


p=\ 

X> 


x 2 


p = 2 
X> 


x 2 


p = 3 
X> 


500-2000 


19 


674 


41.8 


2.20 


23.9 


1.26 


35.1 


1.85 


2000-4000 


17 


901 


84.8 


4.99 


31.2 


1.84 


25.1 


1.48 


4000-6000 


15 


1068 


119.7 


7.98 


33.9 


2.26 


25.7 


1.71 


6000-8000 


13 


697 


109.2 


8.40 


57.2 


4.40 


45.9 


3.53 


8000-12,000 


11 


878 


61.6 


5.60 


22.1 


2.01 


26.4 


2.40 


12,000-20,000 


10 


639 


19.5 


1.95 


19.6 


1.96 


28.2 


2.82 


Entire Sample 


19 


4857 


299.6 


15.77 


111.3 


5.86 


106.9 


5.63 
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10 7 10 8 10 9 10 10 10 11 10 12 



Luminosity (L G ) 

Fig. 1. — Plots of the luminosity function derived by both parametric and non-parametric methods 
for p = 1; the solid curve is parametric, and three non-parametric luminosity functions are shown, 
with different binning in log L. Error bars are suppressed except for the luminosity function with the 
finest binning, which matches the parameterized luminosity function quite closely. The luminosity 
function of SNWZ is plotted for comparison. The histogram at the bottom indicates the luminosity 
distribution of the galaxies in the sample. 
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Luminosity (L G ) 



Fig. 2. — As in Fig. [T], but for p = 2. 
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Fig. 3. — The observed distribution of apparent magnitudes m as a function of log 10 cz. The line 
is the best-fit regression for those galaxies with 2000 < cz < 20, 000 km s _1 . 
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Flux Density (Jy) 

Fig. 4. — Observed (points) and expected (curves) distribution of galaxies in flux density at different 
redshift ranges (labeled in km s _1 ). The solid curves are for p = 1, and the dashed curves for p = 2. 
Topmost box is for the entire sample. Note the virtual indistinguishability of the predictions for the 
two power laws; this shows that the flux density (or equivalently apparent magnitude) distribution 
has little discriminatory power. Predictions are based on the nonparametric luminosity function 
fitted over the range 500-12000 km s _1 , with nu = 25. 
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Flux Density (Jy) 



Fig. 5. — As in Fig. |||, but the luminosity functions are calculated with the SNWZ binning for both 
p = 1 and p = 2. Note that neither curve is a particularly good fit to the data, especially at the 
faint end, although p = 2 does somewhat better. 
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Fig. 6. — The upper panel shows the observed (points) and predicted (curves) distribution of 
galaxies with redshift for p = 1,2,3. The bottom panel shows the ratio of observed to predicted 
counts per bin; the predictions are based on the nonparametric luminosity function. The bar in 
the top panel marks the range of redshifts used in fitting the luminosity function. The points 
in the lower panel show the density field calculated nonparametrically and independently of the 
luminosity function, following Saunders et al. (1990). Solid points are for p = 1, open circles for 
p = 2, and stars are for p = 3. Error bars are shown only for p = 1 to keep the figure from getting 
too cluttered; however, the errors are quite insensitive to the value of p. 
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Fig. 7. — The solid points give the mean log redshift of galaxies in bins of log flux density from the 
entire IRAS 1.2 Jy sample, averaging over galaxies with redshifts between 500 and 20,000 km s . 
The error bars are the error in the mean (that is, the standard deviation divided by the square 
root of the number of points in each bin). The open circles are the mean log redshift for galaxies 
in the Northern Galactic Hemisphere, while the stars are for the Southern Galactic Hemisphere. 
The light solid line gives the expected curve assuming p = 1 and a homogeneous universe, while 
the light dashed line assumes p = 2. The heavy solid and dashed lines give the expected density 
field taking into account the density field found non-parametrically in Fig. n. 
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Fig. 8. — The solid line shows the ratio of the normalized selection function n\(j){z) calculated from 
the interpolated luminosity function described in the Appendix, to that for the analytic luminosity 
function of Eqs. ( |To| ) and (|i~l|). The dashed line is the ratio of ri\(j){z) for the stepwise luminosity 
function to that of the analytic luminosity function. Notice the scalloping in the latter. 



